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Can lattice data for two heavy-light mesons be understood in terms of 
simply two-quark potentials? 
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By comparing lattice data for the two heavy-light meson system (Q 2 q 2 ) with a standard many-body approach 
employing only interquark potentials, it is shown that the use of unmodified two-quark potentials leads to a gross 
overestimate of the binding energy. 
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1. INTRODUCTION 

In this conference the emphasis is on the study 
of single hadrons, with virtually nothing being 
said about systems involving more than three 
quarks. This is unfortunate, since it avoids much 
of particle physics and all of nuclear physics. It 
is, therefore, desirable to attempt to bridge this 
gap. Since the ability to perform accurate lattice 
simulations for systems containing many quarks 
is severely limited - with four quarks being the 
current maximum - the challenge is to see how 
(if) standard many-body techniques can be devel- 
oped for the new situation of interacting quarks 
and gluons. Here we immediately meet a cul- 
tural difference. For particle physicists, the use of 
simply two-quark potentials in multiquark prob- 
lems is not a discussion point, since they know 
(believe) that it is not possible. In spite of this, 
there are many-body physicists, who still believe 
(hope) that this is not so. 

Over the last few years attempts have been 
made to clarify this situation by comparing the 
exact energies of four-quark systems - as calcu- 
lated on a lattice - with standard many-body 
models using only two-quark potentials [1]. The 
four-quark system was studied for two reasons: 

i) Accurate lattice simulations are still possible. 

ii) It is the first step in the description of interact- 
ing hadrons, since the system can be partitioned 
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into two colour singlets. The outcome from those 
pilot comparisons was that - "beyond all reason- 
able doubt" - the resulting binding energies are 
grossly overestimated by the models [Figure [l]] 
and that a four-quark form factor is necessary. 
However, this conclusion was based on lattice 
data for four static quarks in quenched SU(2). In 
the present work most of these approxima- 
tions have been removed by using SU(3) with two 
heavy-light mesons i.e. Q 2 q 2 compared with the 
earlier Q 2 Q 2 . Furthermore, the gauge field is now 
treated in the unquenched approximation. But 
this latter improvement only has a small effect. 
The corresponding many-body model now needs 
some additional assumptions before a compari- 
son with the lattice data is possible. Perhaps the 
most serious is the use of a non-relativistic kinetic 
energy, even though the light quark mass is ap- 
proximately that of the strange quark. However, 
in the binding energy there is a considerable can- 
cellation between the four- and two-quark kinetic 
energies, which then removes much of the kinetic 
energy effect from the binding energy. Also the 
model is developed only for the spin-independent 
contribution to the binding energy. This requires 
an averaging of the lattice data. Both of these as- 
sumptions are not thought to qualitatively effect 
the outcome. The main defect of the calculation 
is that, at present, comparison can only be made 
for small distances between the two static quarks 
- the region where lattice data could suffer from 



2 



0.5 



CO 

CO 



0.4 
0.3 - 
0.2 



> 



CO 

g 

g 
CO 




012345678 

R/a 

Figure 1. The case of four static quarks (Q 4 ) 
in SU(2) placed at the corners of squares of side 
R/a, where a « 0.12 fm. This shows the binding 
energy in lattice units (aB) as a function of R/a. 
The lattice results are for (3 = 2.4 on a 16 3 x 32 
lattice with the dots (crosses) showing the ground 
(excited) state energies - with error bars. The 
two curves are for the model resulting in Eq. ^ 
(with kf = 0) and correspond to the ground (solid 
line) and excited (dashed line) states. This figure 
is similar to the one in 0] 
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Figure 2. Comparison between the spin indepen- 
dent part (Vb) of the Q 2 <f binding energies cal- 
culated on a lattice Q (solid circles - quenched 
approximation with a = 0.170 fm)/ || (solid 
squares - with dynamical fcrmions and a = 0.142 
fm) and the model in the weak coupling limit 
(kf = 0) - solid line. The dashed line shows the 
model prediction for kf = 0.25. These lines are to 
guide the eye. The dynamical fermion data is not 
used in any fit. It is simply included to show that 
it is qualitatively consistent with the quenched 
data but with considerably smaller error bars. 



a lack of rotational invariance. However, in spite 
of these shortcomings, the result is that the use of 
simply two-quark potentials again overestimates 
the lattice binding by upto a factor of three - 
a result that is not expected to be qualitatively 
changed by other models without multi-quark in- 
teractions [Figure || - solid line]. Inclusion of a 
four-quark interaction term can then remove the 
discrepancy [Figure || - dashed line] . 

The reason why the "data-model" comparison 
is good for small R/a in Figure |l| and not in Fig- 
ure H is simply because the Q 2 q 2 energy - even for 
small R - is dominated by potential energy con- 
tributions, where the Qq and qq separations are 
ps 3a. As seen in Figure [l], at such separations 
the model already becomes too attractive in the 
absence of the form factor /. 

2. MODEL 

In Refs. §,§,§ a model was developed for un- 
derstanding the lattice energies of four static 



quarks 

0( r i)Q( r 2)Q(i"3)0( r 4) in terms of two-quark po- 
tentials. This model, in its simpliest form, was 
constructed in terms of the two basis states that 
can be made by partitioning the four quarks into 
two color singlets - namely - 

A = [0(1)0(3)] [0(2)0(4)] 

B = [0(1)0(4)] [Q(2)Q(3)], (1) 

where [...] denotes a color singlet. These two 
states are not orthogonal and have a normalisa- 
tion matrix of the form. 

m = (\ f P). (2) 

In the extreme weak coupling limit the param- 
eter / = 1 and in the strong coupling limit 
/ = 0. However, for intermediate situations it 
is parametrised as 

/(ri,r 2 ,r 3 ,r 4 ) = exp[-b s kfS(rx, r 2 ,r 3 , r 4 )], (3) 



3 



where b s is the string energy, S(ri, r%, r3, is an 
area defined by the positions of the quarks and kf 
is a free parameter. A single value for kf was then 
capable of giving a reasonable understanding for 
100 pieces of data - the ground and first excited 
states of configurations from six different four- 
quark geometries calculated on a 16 3 x 32 lattice. 
In this model the interaction between the quarks 
is expressed as a potential matrix of the form 



V(f) = 



u(13) + u(24) /VC 



AB 



AB 



u(14) +«(23) 



(4) 



where Vab has the form expected in the 
weak coupling limit with the one-gluon-exchange- 
potential 



V = ~^2\iXjv{ij) and v(ij) 



(5) 



i<j 



Away from the weak coupling limit, / is no longer 
unity and in addition v(ij) is taken to be the full 
two quark potential. The energy of the four static 
quarks - a function of the ri - is then given by 
diagonalising 



|V(fe / )-E(4,r i ,fc / )N(fe / )|^ = 0. 



(6) 



This model is easy to generalise to Q 2 q 2 by inte- 
grating over the positions of the two light quarks 
using a suitable variational wavefunction || . This 
results in the need to diagonalise 

|K(*/) + V(fcy) - E{A, R, k f )N(k f )\^ = 0, (7) 

where K(ii, fc/) is the kinetic energy and R 
is the distance between the two static quarks. 
In order to carry out the radial integrals in 
K(R,kf), V(R,k f ) and N(R,kf), the function 
/ in Eq. [| was now reduced to the more symmet- 
ric form 



/ = CX P 



-kfb s 



(8) 



Again kf is a free parameter, which should be ad- 
justed to fit the four-quark lattice energies. How- 
ever, now only the lattice data for R = a is suf- 
ficiently accurate to extract a value of kf — see 
Figure 0. 



The present model is now being developed to 
include: 

1) a semi-relativistic approximation, in which for 
the two light quarks the kinetic energy opera- 
tor in momentum space is simply replaced by 



Pi + m q "V 

2) a spin dependence in the potential between the 
two light quarks. 

This model, although very simple, contains the 
same basic assumptions made in the more elabo- 
rate many-body models that incorporate kinetic 
energy e.g. the Resonating Group Method Q. It 
is, therefore, reasonable that this simplified model 
can to some extent check the validity of its more 
elaborate counterparts. 

3. CONCLUSION 

The comparison between lattice data and a 
model for the Q 2 q 2 system supports (confirms) 
the earlier result - with four static quarks in 
quenched SU(2) - that four-quark energies cannot 
be described simply in terms of two-quark poten- 
tials and that attempts to do so could lead to 
a large overestimate of the binding energy. One 
way of overcoming this problem is to include into 
the model a four-quark form factor - as seen in 
Figure 2. 
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